Asymptotic distributions of periodically driven stochastic systems 
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We study the large-time behaviour of Brownian particles moving through a viscous medium in 
a confined potential, and which are further subjected to position-dependent driving forces that are 
periodic in time. We focus on the case where these driving forces are rapidly oscillating with an 
amplitude that is not necessarily small. We develop a perturbative method for the high-frequency 
regime to find the large-time behaviour of periodically driven stochastic systems. The asymptotic 
distribution of Brownian particles is then determined to second order. To first order, these particles 
are found to execute small-amplitude oscillations around an effective static potential that can have 
interesting forms. 
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I. INTRODUCTION 



When a macroscopic system characterised by a Hamil- 
tonian H(C) is in contact with an environment with tem- 
perature at long enough times it reaches an equilib- 
rium state and is described by the Boltzmann-Gibbs dis- 
tribution P eq (C) ~ exp[— /3H(C)} over the configurations 
of the system. If, however, the system is subjected to 
time-dependent forces of appreciable magnitude, there is 
no analogous general statement that can be made about 
the distribution P(C, t) at large times. A case of partic- 
ular interest arises when these forces vary periodically in 
time. While it is straightforward to see that the large- 
time distribution must be periodic as well, its full form 
is not known in general. It is thus of interest to seek ex- 
plicit answers for particular physical systems subjected 
to periodic driving. 

In this paper, we focus on a paradigmatic system: a 
Brownian particle which feels viscous forces and random 
impulses from the surrounding medium, and is confined 
by a potential well. We ask: What is the effect of a 
further oscillating potential on the state of the particle? 
We are primarily interested in the case where the fluctu- 
ating potential has a nontrivial spatial dependence. We 
analyze the problem mostly in the high-frequency limit 
and find the asymptotic state perturbatively. The re- 
sulting time-averaged asymptotic state is described ef- 
fectively by a distribution of the Boltzmann form with 
an energy function that has three parts: a kinetic energy 
term that depends only on velocity; an effective potential 
that depends only on position coordinates; and a term 
with both velocity and position coordinates. However, 
to the leading order the result is particularly simple, and 
involves only the kinetic term and an effective frequency- 
dependent potential energy whose form can be specified 
exactly. 

The effects of rapidly oscillating periodic forces on 
purely mechanical systems were demonstrated many 
years ago for a driven pendulum [1] and were also studied 
for more general cases [2] . Periodically driven stochastic 



systems too have been studied extensively over the past 
two decades in the context of stochastic resonance (see, 
e.g., Ref. [3], and references therein). However, most 
of these studies were restricted to weak periodic driving 
with position- independent forces [4-7], often only in the 
overdamped regime. Further, it has sometimes been as- 
sumed that the driving frequency is larger than all typical 
frequencies, including that which is associated with the 
noise [5]. Our treatment generalizes that of Refs. [4] and 
[5] in the high-frequency regime, by allowing or arbitrary 
damping; by including position dependence of the driv- 
ing force; and by taking frequencies to be higher than 
those set by the confining potential, yet not necessarily 
higher than those set by the noise. 

Some applications of our results are possible. For in- 
stance, depending on the spatial variation of the fluctu- 
ating forces, the effective potential may have more than 
one local minimum even if the original potential has only 
one. Under such circumstances, an assembly of Brow- 
nian particles would tend to segregate in two separate 
collections. Moreover, the fact that the effective poten- 
tial depends on physical properties, such as the mass of 
the particles, can be exploited to promote segregation of 
two sets of Brownian particles that differ from each other 
in mass or some other physical attribute. 

The layout of the paper is as follows. In the following 
section, we define the problem and discuss the different 
time scales involved and their interplay, and discuss var- 
ious regimes qualitatively. In Sec. Ill, we develop the 
necessary formalism required to address rapid periodic 
drive and arrive at a perturbative scheme. In Sec. IV, 
we use this scheme to determine the asymptotic distri- 
bution. In Sec. V, we briefly discuss the effect of slow 
periodic driving so as to compare with that of rapid driv- 
ing. Finally, in Sec. VI, we conclude with a discussion 
of possible directions in which our results may be gener- 
alized, possible applications (particle segregation, parti- 
cle sifting), and an example (the simple pendulum) that 
demonstrates the significance of the new effects found. 
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II. PERIODICALLY DRIVEN BROWNIAN 
PARTICLE 

We shall consider a one-dimensional Brownian particle 
in a potential well, subjected to a periodic force along 
with a damping force and random noise. The equation 
of motion of the driven Brownian particle moving in a 
viscous environment is 

d 

mx = -jx - g^U{x) + F(x,t) + T](t), (1) 

where m is the mass, 7 is the coefficient of viscosity, 
U(x) is a static confining potential, F(x, t) is the peri- 
odic driving force with a period T, F(x, t) — F(x, t + T), 
and r](t) is a Gaussian random noise with (T](t)) ri = and 
{v(t)v(t'))r, = 27/3- 1 5(t - t'), where /3" 1 is the temper- 
ature of the surrounding heat bath. 

The probability distribution of the Brownian particle is 
defined as P(x, v, t) = (5(x — x v {t))5{v — x^t)))^, where 
x v (t) and x v {t) are the position and the velocity at time 
t for a particular history of {ri(t)} over a time t. The 
time evolution of P(x, v, t) is described by the follow- 
ing Fokker-Planck (FP) equation (also referred to as the 
Kramers equation), 



dP _ d_ d_ 

dt dx dv 



— {-1v-U'(x)+F(x,t)}P 
m 



+ 



7 d 2 P 
(3m 2 dv 2 



(2) 



The driving force F(x,t) that is oscillating with a fre- 
quency u) — 2n/T is chosen to be 



F(x,t) = f(x) cos(ujt) + g(x) sin(wf) 



(3) 



in the domain L\ < x < L2, where the amplitude func- 
tions f(x) and g(x) vanish at the boundaries, x = L\ 
and L 2 , and outside the domain. This choice of F(x,t) 
is made for convenience; choosing a more general peri- 
odic function will not hinder our analysis. The gener- 
alization to higher dimensions and to many interacting 
particles is also straightforward. Our aim is to find the 
large-time distribution that we denote by P OQ (x,v,t) = 
finite P(x,v,t). 

In the absence of a driving force, all solutions P(x, v, t) 
of the FP equation, corresponding to various arbitrary 
initial distributions, tend to a unique distribution af- 
ter a long enough time [8]. This distribution P 00 (x,v,t) 
for the Brownian particles takes the equilibrium canon- 
ical form P eq (x,v) = (I/Zq) exp {— (3 [^mv 2 + U(x)] }. 
When a periodic driving force is present, Poo(x,v,t) ap- 
proaches a periodic function of time which is unique up 
to a phase [9,3]. In brief, the argument for the peri- 
odicity goes as follows. When the FP operator is peri- 
odic, C(t) = C(t + T), the solution to the FP equation 
[dt — C(t)]P(x, v, t) = can be expanded in terms of the 
Floquet-type functions p^{x,v,t)ex^>(—^,t). The func- 
tions p^(x,v,t) = Pfi(x,v,t + T) are periodic and are the 



right eigenfunctions of dt — C(t) with eigenvalues /i. It is 
known that for an iV-dimcnsional FP equation, the real 
parts of these eigenvalues Re(/i) are positive semidcfinite 
and hence in the large-time limit, for typical confining 
potentials, only p (x,v,t) survives. We are interested in 
finding this large-time distribution P OD {x 1 v, t) for a given 
F(x, t). Since no analytic solution of the FP equation is 
known for an arbitrary time periodic F(x, t), even when it 
involves only the fundamental frequency, we shall restrict 
our attention to certain regimes of the driving frequency 
while solving for Poo{x, v, t). 

In the absence of the driving force there are two impor- 
tant time scales in the system; one, t v = m/7, is intro- 
duced by the viscous medium and the other, t w — 2tt/ljo, 
is related typically to the curvature at the bottom of 
the potential well, loq = ^U"(x m i n )/m. The veloc- 
ity variable equilibriates in a time scale set by t v and 
hence for larger times it gets described by a station- 
ary distribution. Thus, in a highly viscous medium, 
where t v <C t w holds, after a time t 3> r v the distribu- 
tion can be written as P(x,v,t) = P eq (v)P(x,t), where 
Peq(v ) ~ exp [— P^mv 2 ] is the canonical distribution for 
velocity and P(x,t) = J dvP(x,v,t) is the marginal dis- 
tribution involving only position. This marginal distri- 
bution satisfies an FP equation (often called the Smolu- 
chowski equation) obtained by dropping the inertial term 
in Eq.(I), as this term becomes insignificant in compari- 
son with the viscous term once the time exceeds t v . 

If, however, there is a driving force with a high enough 
frequency then this reduction to the Smoluchowski limit 
does not take place. The oscillating driving force in- 
troduces one other important time scale associated with 
the timeperiod T whose existence restricts the domain of 
validity of the Smoluchowski equation even in the high- 
friction limit. It is useful to demarcate different regimes 
of the driving frequency in this limit: (a) t v <§C t w < T 
and (b) t v < T <ti t w or T < t v <ti t w . In case (a), the 
velocity decouples from the position and one can still use 
the Smoluchowski equation to obtain the large-time dis- 
tribution, while in case (b), the driving force has a time 
scale comparable to that of the velocity relaxation time 
and hence it is necessary to retain the Kramers equation. 

If the time scale of measurement exceeds the time 
period T, then the relevant quantity is t he large-tim e 
distribution averaged over a time period, P OQ (x,v,t) = 
(l/T) J Q T Pocix, v, t). Hence we shall also determine 

Poo(.T, v,f). 



III. RAPIDLY OSCILLATING FORCES: 
FORMALISM 



In a mechanical system, in the absence of a viscous 
force and random noise, it is known that for w » wj, 
the particle executes small amplitude oscillations of fre- 
quency iv about a smooth mean path [2]. This motion 
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can be described by separating it into slow and fast vari- 
ables; the fast variable decouples from the slow one while 
the slow variable sees a static effective potential, modi- 
fied due to the oscillations of the fast variable. In this 
section, we develop the necessary formalism that accom- 
modates an analogy with this separation into slow and 
fast variables and enables us to solve the problem when 
viscous and random forces are included. 



where f v {X) = -dU(X)/dX. 

Now choose £(A, r) such that it is a solution of the 
following equation 



m£ = -rt + F(X,T) . 



(7) 



For F(X, t) as given in Eq. (3), the solution to the above 
equation is 



A. Transformation of FP equation 

In this subsection, we transform the FP equation under 
a specific coordinate transformation which enables it to 
be solved perturbatively. For the perturbative treatment 
to be valid a sufficient condition, though not necessary, 
is that u) be large, while no assumption is made about 
the amplitude of the driving force in comparison with 
the static potential. 

We make the coordinate transformation {x,v,t} — > 
{X, V, t} under which the distribution is made to be- 
have like a scalar function: P(x,v,t) — > P(X,V,t) = 
P(x,v,t). The old and new coordinates are related as 
follows, 

x = X + £{X, r), v = V + £-£(X, r), * = r. (4) 

The explicit form of £(A, t) will be specified later. Note 
that the volume element dxdv = dXdV J(X, V, r), where 
J(X, V, r) = (1 + dZ/dX) and hence, if P is the probabil- 
ity density in (x, v) space then P J becomes the the prob- 
ability density in (X, V) space. Also, under the above 
coordinate transformation the derivatives transform as 
follows: 



d 



1 



d 



i' 



d 



dx l + f dX l + f dV 
d_ _ d_ 

d_ _ d_ _J d_ ( jj' 

di~!fr~ l + f dX + ( 1 + f 




(5) 



where the dots and primes on £ stand for derivatives with 
respect to r and X, respectively. Making use of Eqs. (4) 
and (5) in Eq.(2), we get 



dP 



dx K ' av 



i 



rri 



{--rV + Mx + s) 



+ v 
i 



F(X + £,t)-F(X,t)}P 
f dP g dP 



7 d 2 P 



+ 



m 



i + ?dx i + ?dv 
dP 



m£ + rt-F(X,T) 



(6) 



-1 



/(A) + -^g(X) ) cos(wt) 



+ [ g (X)--Lf(X))sm(u;T) 
mui 



(8) 



Since £ is small for large values of (to 2 +7 2 /m 2 ), we may 
expand Eq.(6) perturbatively in £. The reason for choos- 
ing the above dynamics (Eq.(7)) for £ is that the last 
term in Eq.(6) becomes zero and further the explicitly 
(time)r-dependent term becomes small if £ is small, thus 
making the equation amenable to perturbative analysis. 

Note that when the driving force is x independent 
Eq.(6) reduces to the usual FP equation with the static 
force /[/(A) being replaced by fu{X + £). 



B. The perturbative scheme 

We now find the large-time solution of Eq.(6) pertur- 
batively in powers of £. Upon substituting £ as given in 
Eq.(8), we see that Eq.(6) takes the form 



a7 



P(X,V,t) = [C + AC]P(X,V,t) 

oo 

= ^[£ (n) + A£W]P(X,V,t), 



(9) 



n=0 



where C and A£ are static and time-dependent opera- 
tors, respectively. The superscript on the operators 
and A£(") indicates that they are of 0(£ n ); the explicit 
forms of the first few operators are 



£ (0) 



aX m oV pm z aV z 



I 







-rf^r) 8V' 



A£(°) = 0, 

A£« = - 
m 



ZU"(x)-i;f , (x,t) + zf>(x,t) 



dV 



A£< 2 > = J- 



(e-e)u'"(x)-eF"(x,r) 



1 

2m 

l? 5 >dx ^ av 



d_ 

dV 



(10) 
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where the bar over the terms indicates an average over 
a time period. The perturbative asymptotic solution of 
Eq.(9) can be formally written as follows 

P x (X, V, t) = Qoo (X, V, t) + -^—^ AC Poo (X, V, t) , 

(11) 

where Q oa (X, V, r) is the right eigenfunction of C with 
eigenvalue zero. Since the asymptotic distribution is pe- 
riodic and the nonzero eigenvalues of C have a nonvan- 
ishing real part, it follows from Eq.(9) that A£Poo has 
no overlap with the eigenfunction of d T — C with zero 
eigenvalue. Hence the operation of the inverse of d T — C 
is well defined as it acts only on the space of functions 
orthogonal to that eigenfunction. 

It is not possible to obtain an explicit form for the in- 
verse of d T — C for arbitrary U(X) and F(X, r), in which 
case numerical or variational methods might be adopted 
to determine the eigenfunctions and eigenvalues of C or 
Co and to obtain (d T — £) _1 . If for some form of the po- 
tential U(X), we are able to find all the eigenfunctions of, 
say, Co then this inverse can be conveniently expanded 
to any order in £ as follows 



w oc 



n=0 



1 



d T - C 



(C - Co) 



1 



9 T - C Q ' 



(12) 



But a more suitable series expansion can be given for this 
operator when either u> or 7/m is large. If to is also much 
larger than 7/m, then this suitable expansion is 



1 



= E 



n=0 



C 



1 



(13) 



In the case when 7/m is comparable to w, we can write 
the inverse in terms of an operator Cy containing terms 
of 0(j/m) only and which is given as 

Cv = + ^L_S^ = ^ mV2 f--«M e*' 

V m / 



m dV 13m 2 dV 2 



where the operators a and follow the commutation 
relation [a, a^] = 1 and are defined as 




(15) 

So in this case the convenient expansion of the inverse 
operator is 



1 



d T -C 



E 

n=0 



e * 



1 



d T - C\ 



1 



d T - C\ 



The idea of writing this operator in terms of a and 
is that its action on h(V) exp(— f3mV 2 /2), where h(V) 
is some polynomial of V, can be determined more easily 
since it involves the action of a specific function of a^a 
on a series made of eigenfunctions of a) a. 

The calculational scheme is thus reduced to the follow- 
ing: Find the right eigenfunction Q^ of C with eigen- 
value zero to the desired order in £. Then use one of the 
truncated series (12), (13), or (16), such that the trunca- 
tion is consistent with the chosen order, keeping in mind 
that £ depends on ui. Next substitute this truncated in- 
verse operator in Eq.(ll) and from it extract P OQ (X, V, r) 
order by order. 



IV. ASYMPTOTIC DISTRIBUTION: 
PERTURBATIVE ANALYSIS 



A. First-order perturbation 



We now find the asymptotic distribution P oa (x,v,t), 
to first order in £, with the condition that it vanishes 
at the boundary of (x, v) space. We will also show that 
the time average of this distribution P 00 (x,v,t) is the 
canonical equilibrium distribution with U(x) replaced by 
U eff (x) = U(x) + UW(x), where U {1 \x) will be explic- 
itly evaluated. 

To this end we need to determine Qoo (X, V, t) = 
Q^oo {X, V, t) to the same order from the equation 
CQoo{X, V, t) = 0. As before, the superscript on Q^ 
indicates the corresponding order in £. The zcroth order 



Qoo is the solution to the equation C^QrJ — which 
yields the equilibrium distribution 



(14) P£\X,V,t) = QW(X,V,t) 



= ^- exp ( - >' 



\mV 2 + U(X) 



(17) 



The first-order Pi^ is obtained from Eqs.(ll) and (16) 



1 e^ 2 A£« PW, 



d T + —a^a 

1 1 , m 



(18) 



x „ where Qoo is the solution to the equation C^QiJ + 

(- ^y—f- - £v)) ^^Qi? = 0. This solution is straightforward to deter- 

n=o \ T ^ m a a mine since C^Q^} has the same form as the right hand 

side of the identity: £<°) [h(X)Q l ) ] = -Vh'{X)Q { £ for 
any arbitrary function h{X). We get 



(16) 
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Qg\X,V,T) = -0UW(X)pW 
d 



dX 



uW(X) = -Z{X,t)F'(X,t) 



(19) 



In Eq.(18) we have used the inverse operator (16) trun- 
cated right after the first term which amounts to neglect- 
ing 0(£/u)) terms. To this approximation, it is then suf- 
ficient to keep only the term proportional to £ in AC^\ 
which is rewritten as follows 



„ 9 tt + a' •., x. t 

+ ^Wr~R £(« + <*V 

oX [5m 



e 4 



(20) 



Hence Eq.(18) upon neglecting terms of 0(£,/uu) reduces 
to 



P£\X,V,t) = -(3UW(X)Pg> 



1 



-^mV 2 }_ p -()U(X) 

Z 



i{a + a t )a t 

p[u^(X)+K^(V;0\P^(X,V,r), 

(21) 



where 



/3m \ m 



2H 



27;, 



m 



£ mV 



(22) 



which simplifies to K^\V;^) s=s £mV 2 when uj 3> 7/m. 

We now get the asymptotic distribution from Eqs. (4), 
(17), and (21), 



Poo(x,v,t) = — cxp( -/? 



-m» 2 + C/(a;) + [/W(x) 



+tf( 1 )( U ;£) -£[/'(*)-£ 



mu 



(23) 



Unlike in the distributions for static potentials, here the 
velocity v gets coupled to the position x through £. The 
averaged large-time distribution is given as 



1 



P<x,(x,V,t) = — cxp( -f3 



\mv 2 + U{x) + U (1) {x) 



(24) 



where the explicit expression for (x) is obtained upon 
substituting £ in Eq.(19) and then integrating: 



U (1) (x) = 



1 



4to(w 2 + Xj) 

27 



f{x)+g\x) 



+ [ X dy(g(y)f(y)-f(y)g'(y) 
muj J \ 



• (25) 



Thus the time-averaged large-time behavior of the Brow- 
nian particles is described by the canonical distribu- 
tion at a temperature /3" 1 with an effective potential 
U e ff = U + that depends on the frequency and 
space dependence of the driving force in addition to the 
properties of the particle. Note that a nontrivial contri- 
bution to the effective potential arises to this order only 
if / or g are space dependent, and that it can be tuned 
by varying /, g, or u). 

The additional term is the average energy asso- 
ciated with the rapid motion. This can be seen easily 
upon substituting for F in terms of £ using Eq. (7) and 
rewriting as follows: 

uW(x) = ~ f dy Uy,t)F'{ y ,t)] 



= \m [i(x, t)] 2 - f dy (-7^, t) ^(y, i)) . (26) 

Thus is sum of the average kinetic energy and the 
work done against the damping on the fast variable £. 

A first-order perturbation treatment is justified, pro- 
vided £( 2 ) and AC^ are negligible when compared to 
and AC^\ respectively. The criterion for this is 
\t;df/dx\ » |£/"'£ 2 |. For instance, suppose that the 
length scales over which U(x) and F(x, t) vary are com- 
parable and au>o is a typical frequency associated with 
anharmonic terms. Then the above criterion reduces to 
(iv 2 + 7 2 /m 2 ) 3> a 2 uj 2 , which is consistent with £ being 
small. 

It should be remarked that the restriction of period- 
icity of F(x, t) to the fundamental frequency u> is not 
essential; one can include the higher harmonics as well. 
Also by considering generalizations to higher dimensions 
or to many interacting particles one does not encounter 
any additional computational difficulties in the evalua- 
tion of the distribution. The only additional assumption 
needed to write down U e jf is that f(x) and g(x) are curl 
free. 



B. Second-order perturbation 

We now calculate the second-order corrections to the 
asymptotic distribution. Further, in this subsection, we 
restrict our treatment to frequencies that satisfy the con- 
dition oj 3> j/m, in which case the first-order term P^' 
is as given in Eq.(21) with K^(V;C) = £mV 2 . The 
second-order term of Eq.(ll) is 



= Q { 3 + 



[A^WpW + AC^PV + A£( 2 )PW] 



A£ (l)p(0) 



^AC^PW 



+±ac^ p^+o(-e)- 



UJ 



(27) 
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Though the second term in the above expression is of 
0(£) it has been included here because when calculated 
earlier, to first order in £, the terms of 0(£/uj) were ne- 
glected. This term is written within pipes to indicate 
that only terms of 0(£/u) and 0(^/u 2 ) have to be re- 
tained and not the 0(£/cu°) term that has already been 

included in P^ . 

We now evaluate the terms on the right-hand side of 
Eq.(27). The first term Qio is the solution to the equa- 
tion £(°)Q$ + £«q£ } + C^Q ( £ = 0. This, as in the 
case of , is straightforward to determine and we get 

Q&\X,V,t) = -PUi(X)PW + 1 (u^{x)) 2 P£\ 



The second term is 
1 



1 



(aru'(x) + -(tfu'"(x). (28) 



A£ (D p{0) 



= -/3/C(X,F;0 PW, 



(29) 



where 



K{X, V; Pi 0) = P(X, V; - ±£W H(X, V; £), 
H(X, V; = (tU" + t?U' - \(iF> + VP<® 
+£M (mF^'P^) . (30) 
The third and the fourth terms are 
I. A£ (D pW = -pmV^'f W]P£ ] 



+/3 2 [[/( 1 W 2 £' + L(mV*tr]P£ ) , 



1_ AC (2) p(0) = ^ m y2 [( ^ )2 _ (e)2] p(0) 

-/3[W 2 (X) + C(X,t/)]Pi o ), 



(31) 



where (W 2 )' = tt'U" - 3£' £' U' and C(X, V) is an ar- 
bitrary r-independent function. These arbitrary and r- 
independent terms are included since the action of d~ x 
allows this ambiguity. This ambiguity is removed from 
the condition obtained by substituting P^ in Eq. (9) and 

~ (2) 

averaging over a period. Thus PV should satisfy 



c<npg ) + A£WP£ ) =o, 



(32) 



which implies, after some calculation, that C is the solu- 
tion of the equation 



£(°)[CPW] + 4m^"V 3 P^ = 0. 



In the high- viscous limit, j/m 3> luq, this function 
C(X,V) w (4m/37)FF(m^ 3 + 6V/0). Using Equa- 
tions (27)-(31) and Eq.(21) we get the asymptotic dis- 
tribution to second order 



Pco(x,v,t) = |exp(-/3 \mV 2 + U(X) + U^(X) 



+U^ 2 \X) + R{X,V,t) 



(34) 



where R(X, V, t) 
\mV 2 [(S,'(X,t)) 2 



_^ K^(V;0 + £(X,V;0 + 

(£W)) 2 ] +C(X,V); U^(X) = 
U\{X) +U 2 (X); and, to second order, X = x — £(x, t) — 
£(x,t)£'(x,t) and V = v - £{x,t) - £(x,t)£'(x,t). It is 
clear from the above equation that P 00 (x,v,t) will con- 
tain terms with both v and x dependence in addition to 
purely x-dependent and w-dependent terms. 

The explicit form of U^ 2 \X) = U^ 2 \x) which con- 
tributes at second order to the effective static potential 



is 



4m 2 to 2 (to 2 + X) 



f(xf+g(x) 2 )U"(x) 



dy[(.r(y)) 2 + ( 9 '(y)) 2 )U'(y) 



(35) 



(33) 



We notice here that a nontrivial contribution arises even 
for x-independent driving provided U (x) is anharmonic. 
Also note that this second-order correction to the effec- 
tive potential depends on U (x) while this is not the case 
at the first order. 

We recover the results of Dcvorct et. al. [4] and Jung 
[5] when f(x) and g(x) are independent of the posi- 
tion coordinate x. In the former reference the iner- 
tial term was considered while in the latter it was not, 
but their analysis of the high-frequency limit is tanta- 
mount to assuming that the driving frequency is larger 
than all typical frequencies of the system including that 
of noise. This assumption restricts the validity of the 
answer, and clearly does not hold for white noise. In 
fact, this assumption would lead to the absence of the 
term J [(f) 2 + (g') 2 ]U' in the effective potential for re- 
dependent driving. That the error, when this assumption 
is made, shows up at 0(£ 2 ) is also evident from the pres- 
ence of the term V£' 2 d/dX in operator. 

We reiterate that there is a nontrivial contribution to 
the effective potential from the position-dependent driv- 
ing forces. First, it shows up at first order itself and, sec- 
ond, it shows up at second order even for harmonic poten- 
tials; both of which are absent for position-independent 
driving forces. 
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V. SLOWLY OSCILLATING FORCES 

In this section, we will consider the opposite extreme, 
namely, oj small compared to both 7/m and u>q. The aim 
is to compare with the asymptotic behavior we found 
in the previous section under rapid driving. Under slow 
driving, the Brownian particle sees an unchanging poten- 
tial within its relaxation time and so the Boltzmann-like 
distribution corresponding to the instantaneous potential 
is a good zeroth-order starting point for perturbation. 

The large-time distribution Poo(x, v, t) will acquire the 
same periodicity as that of the driving force. Hence 
the left-hand side of Eq.(2) is of O{to) while the right- 
hand side has two terms, one of 0(-f/m) and the other of 
0(loq)- Thus P 00 (a;,f,t) can be evaluated perturbatively 
in (mw/7) and (uj/ojo). To the leading (zeroth) order, 
this distribution is 



exp{-/3[imu 2 + U {x) 
+U f (x) cos(ui) + U g (x) sm(ut))]}, (36) 



where U f (x) = - J x dyf(y) and U g (x) = - J x dyg(y) ; 
Z{t) is the instantaneous normalisation constant deter- 
mined by the normalization condition: J dxP^ (x, v,t) = 
1 at any time t, and evidently satisfies the periodicity 
condition Z(t + T) = Z(t). 

We now estimate the averaged large-time distribution 
P oa (x,v,t) at low and high temperatures separately. To 
determine the average distribution at high temperature, 
it is convenient to rewrite this equation as follows. 



n#0 



Jn (PV(X)) 

Jo ((3V(x)) 



(37) 



where X„(a) is the modified Bessel function and V(x) = 
y/[U f (x)] 2 + [U g {x)\ 2 . The ratio l n (a)/l (a) lies be- 
tween and 1 for < a < 00 , and decreases very rapidly 
for small a; l n (a)/Xo(a) ~ a n /2 n n\. Thus it suffices 
to keep a few terms, and to leading order we obtain 
Poo{x,v,t) = (l/Z) exp{-/3[±mv 2 + U eff {x)]}, where 

U eff (x) = U(x) -- p \n{l {^[U f (x)] 2 + [U g {x)Y)}. 

(38) 

Thus the Brownian particles, when observed over a time 
of 0(T), get described by the canonical distribution at a 
temperature with an effective potential that depends 
on the temperature. 

At low temperature the saddle-point approxi- 
mation can be used to evaluate Poo(x,v,t). To 
the leading order this will be Poo( 



P eq (v) [l/N(t)\Y^?5{x-xZ n {t)), where 



W 

rain \ 



,{x,v,t) 
.(0 



are the minima at time t of the function U(x) + 
Uf(x) cos(wt) + Ug (x) sm(cot) and N(t) is the number 
of minima at that time. 

It might appear that U f(x) and U g (x), obtained by in- 
tegrating f(x) and g(x), respectively, have different con- 
stants of integration in regions where the driving forces 
are present from those where they are absent. How- 
ever, they get fixed by the condition that P(x,v,t) and 
dP(x, v, t)/dx are single- valued functions of x. 



VI. DISCUSSION 

In this section, we review the results of this paper for 
a general position-dependent periodic driving force with 
a focus on the high frequency regime. We then discuss 
effective potentials for specific forms of the driving func- 
tions, and use these results to point out some possible 
applications. The nature of these effective potentials can 
be drastically different from that of the original poten- 
tials, and we illustrate this fact using the example of a 
simple pendulum. 

The principal results of this paper concern the form of 
the asymptotic distribution of a Brownian particle sub- 
jected to a driving force, which is periodic in time but 
is an arbitrary function of position. In the limit of high- 
frequency driving, the particle makes small, rapid excur- 
sions around a smooth path along which the motion is rel- 
atively slow. This affords the possibility of a systematic 
perturbative treatment in powers of the excursion am- 
plitude. The result for the distribution averaged over a 
cycle is described in terms of an effective potential whose 
form we derived. Interestingly, the leading contribution 
U^ 1 ' , which is second order in the amplitude of the ap- 
plied driving force, is present only if the driving force is 
position dependent. This can be interpreted as the 
average energy of the excursion variable. In the next- 
order contribution to the effective potential as well, this 
position dependence of the driving force has an interest- 
ing effect even for purely harmonic confining potentials. 

Central to our discussion of rapid periodic driving is 
the separation into slow and fast variables. It should be 
noted that this demarcation of slow and fast is based on 
whether or not the variable varies considerably over a 
time period. This is different from the distinction made 
in discussions of fast and slow variables that are so cate- 
gorized according to whether they relax in small or large 
times. In this latter case, many methods have been de- 
veloped to eliminate the fast variables and obtain an ef- 
fective dynamics for the relaxation of the slow variables 
[8,10]. 

The effective potential that we find can be qualitatively 
different in the low- and high-frequency regimes. When 
cj is small, the leading additional potential felt by the 
particles, AU(x) = U e ff(x) — U(x), is always nonposi- 
tive for any choice of f(x) and g(x). On the other hand 
for large u, AU(x) can be either positive or negative in 
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general, though if one of f(x) or g{x) is identically zero 
then it is always non-negative. Also to the leading order 
when w is small, AU(x) depends on temperature but not 
on co, whereas in large to limit it only depends on u and 
not on temperature. 

Certain choices of the driving force can lead to in- 
teresting outcomes. Suppose we choose U(x) to be a 
confining potential which is a monotonically increasing 
function of |x| while f(x) is a monotonically decreas- 
ing function of |x| which vanishes at \x\ — L, and g(x) 
vanishes identically: for example, U(x) = ^mufix 2 and 
f(x) = /o sin(7nc/L) for — L < x < L and zcroothcrwise. 
When co is small the effective potential continues to have 
a single minimum at the origin. But when u> is large, it 
develops two additional minima; one close to L, between 
and L and the other close to — L, between — L and 0. 
Thus, a system of Brownian particles would cluster in a 
region near the origin for low driving frequency, while at 
high frequency these particles would segregate into two 
clusters which are separated by a distance of O(L). 

The parameters specifying the particles also enter the 
effective potential. Hence the minima of the effective 
potential as seen by different species of particles are dif- 
ferent. One can make use of this fact, for instance, to 
separate different species of Brownian particles in a sit- 
uation where they are initially mixed, by driving them 
with a space-dependent periodic force. 

To get an idea of the magnitude of the qualitative 
change the effective potential introduces it is instructive 
to examine an example. Consider a rigid pendulum: a 
massless rod of length I with a bob of mass m attached 
to it at the end and oscillating in a gravitational field g. 
The potential has extrema at 9 = and 9 = tt which 
are respectively stable and unstable points, where 9 is 
the angle the rod makes with the negative y axis (along 
the direction of gravity). Now oscillate the point of sus- 
pension along the y axis with a frequency to and am- 
plitude a. The angle 9 then evolves according to the 
equation 19 + g sin 9 = au 2 sin 9 cos(wt). The effective 
potential V e ff(9) has extrema at 6 = 0, tt, and 9±, where 
9± = tt ± cos -1 A with A = 2gl/a 2 uj 2 . The stability of 
these points is as follows: 9 = is stable; 9 = 9± ex- 
ist only if A < 1 and when they exist they are unstable; 
9 = tt is stable when 9± exist and is unstable otherwise. 
In a nutshell, when A > 1 the pendulum shows no qual- 
itative change in its behavior upon oscillating the point 
of suspension whereas when A < 1 then 9 = tt also be- 
comes a stable point and hence the pendulum can make 
small oscillations about this point too. This dramatic 
change in the behavior of the pendulum was experimen- 
tally demonstrated by Kapitza [1]. 

Interestingly, if we oscillate the point of suspension 
along the x axis instead of the y axis, the pendulum will 
exhibit a different behaviour. In this case, the effective 
potential has extrema at 6 = 0, tt, and ± cos^ 1 A if A < 1; 



the points ± cos -1 A are stable when they exist; 9 = is 
unstable if A < 1 and stable otherwise; 9 = tt is always 
unstable. In other words, the pendulum now oscillates 
about points that do not lie on the y axis. The nature 
of the extrema of the effective potential does not change 
in the presence of the viscous term and noise except that 
the value of A will now be different (with to 2 replaced by 
^/(u; 2 + ( 1 /m) 2 )). 

To conclude, we have developed a perturbative cal- 
culational scheme to study the asymptotic behavior of 
Brownian particles under the influence of rapidly oscil- 
lating forces. When these forces are position dependent, 
nontrivial effects are seen in the large-time behavior. The 
formalism developed here can be generalized straighfor- 
wardly to interacting Brownian particles in any dimen- 
sions. It can also be extended to study the behavior of 
fluctuating fields when subjected to periodic driving. 
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